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<D ■ 
^ ■ Abstract 



We study finite temperature J\f = 4 Super Yang-Mills (and more general gauge theories 
realized on intersecting D3-D7 branes) in the presence of dynamical massless fundamental 
matter fields at finite baryon charge density. We construct the holographic dual charged 
black hole solutions at first order in the flavor backreaction but exact in the charge density. 
The thermodynamical properties of the dual gauge theories coincide with the ones found in the 
usual charged D7-probe limit and the system turns out to be thermodynamically stable. By 
analyzing the higher order correction in the flavor backreaction, we provide a novel argument 
for the un-reliability of the charged probe approximation (and the present solution) in the 
extremality limit, i.e. at zero temperature. 

We then consider scalar mesonic-like bound states, whose spectrum is dual to that of 
linearized fluctuations of D7-brane worldvolume fields around our gravity backgrounds. In 
particular we focus on a scalar field saturating the Breitenlohner-Freedman bound in the 
flavorless limit, and coupled to fields dual to irrelevant operators. By looking at quasinormal 
modes of this scalar, we find no signals of instabilities in the regime of validity of the solutions. 
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1 Introduction 

The study of the Quantum Chromodynamics (QCD) phase diagram requires a full compre- 
hension of finite temperature (T) and baryon chemical potential (//) regimes which are often 
beyond the reach of current theoretical methods. Standard perturbative treatments of QCD 
only provide access to certain asymptotic corners of the (T, //) plane. Lattice simulations are 
limited by the so-called sign problem: a finite value of /i leads to unreliable Monte Carlo 
simulations (though recent advances, e.g. with analytic continuations to imaginary potential, 
are partially improving the situation). Moreover, since they are based on a reformulation of 
QCD on discretized Euclidean spaces, lattice methods are not well suited to study real-time 
processes. 

A challenging phase of QCD is the Quark-Gluon Plasma (QGP), in which the nuclear 
matter dissociates into its fundamental components: quarks and gluons. This state of matter 
is created in laboratories like RHIC and LHC by collisions of heavy nuclei. The dynamics of the 
QGP is described by strong-coupling interactions and involve a finite (though typically small 
if compared with the plasma temperature) quark density. Moreover, most of its interesting 
properties reside in its real-time dynamics (hydrodynamic transport coefficients, quenching 
of probes moving through it, etc). Standard theoretical methods are therefore not exhaustive 
to comprehend its features. 



Another particularly interesting part of the QCD phase diagram involves a "large" chemi- 
cal potential (at least comparable to the zero charge crossover temperature) at zero and finite 
temperature. This is the region where a critical point, if it exists, should reveal itself. Due 
to its strongly coupled and large \x nature, this region of the phase diagram is beyond any 
current first-principle theoretical exploration. 

A tool that has been used long and wide in the last years to understand the dynamics of 
strongly coupled quantum field theories is the gauge/gravity (holographic) correspondence. 
This allows to map classes of strongly coupled quantum field theories into weakly coupled 
theories of gravity in at least one dimension more. The extra (radial) direction is the geometric 
counterpart of the Renormalization Group (RG) scale of the dual field theory. 

The best understood realization of the correspondence, known as AdS/CFT, relates grav- 
ity theories on asymptotically Anti-de-Sitter (aAdS) backgrounds in d + 1 spacetime dimen- 
sions to conformal field theories in d-dimensional Minkowski spacetime. Finite temperature 
and chemical potential phases in these theories are then mapped into dual aAdS charged 
black hole backgrounds. 

Although there are several examples in which the correspondence between a quantum field 
theory and a dual theory of gravity is explicitly realized, this is unfortunately not the case for 
QCD: the holographic dual of Quantum Chromodynamics remains out of reach. Maximally 
super symmetric (A/" = 4) SU(N C ) Yang-Mills (SYM) theory in four dimensions, instead, has 
a well known holographic dual, and most of its non-perturbative properties can be exactly 
deduced using the correspondence. It is a conformal field theory with a bosonic content of six 
scalars and one vector in the adjoint representation; the scalar manifold has 50(6) symmetry. 
Despite the enormous differences among this theory and QCD, in the last years we have 
learned that it provides sensible benchmarks on some equilibrium and dynamical features of 
phases where QCD could be roughly approximated by a (strongly coupled) conformal model. 
This is the case of the deconfined QGP phase, where e.g. the trace anomaly results to be quite 
suppressed for not too large values of the temperature above the confinement-deconfinement 
crossover. On the other hand, turning on temperature badly breaks supersymmetry, thus 
allowing a phenomenologically unpleasant feature of the N = 4 SYM model to be cast away. 
See [lj for a recent review on all these issues. 

A QCD-inspired refinement of the above toy model requires at least the inclusion of quarks, 
i.e. fermionic matter fields transforming in the fundamental representation of the gauge group. 
Actually, in order to preserve some supersymmetry in the T = vacuum (so to benefit from 
a quantitative control on the related gravity dual) one has to introduce flavor multiplets, 
containing scalar (s-quarks) fields as well. The inclusion of these multiplets generically breaks 
part of the global symmetries of the original model. The case we will focus on in this paper is 
that of massless fundamental hypermultiplets which, in the N = 4 SYM case, are introduced 
in such a way to preserve an 50(4) x U(1)r C 50(6) symmetry subgroup, where U(1)r is 
an R-symmetry. 

The addition of fundamental degrees of freedom also breaks conformal invariance, though 
in an unpleasant way. The beta function for the 't Hooft coupling A = gyM^c, results to be 
positive (thus signaling the occurrence of a UV Landau pole) and proportional to the number 
of flavored species, /3(A) ~ (Nf/N c )X 2 , as opposed to the negative beta function of QCD. 
Just as in the case of QED, another famous theory exhibiting a UV Landau pole, a sensible 
treatment of the model can be obtained only by focusing on the IR physics. 

The present work will actually focus on M = 4 SYM theory (and its "quiver" gener- 
alizations thereof, see below) coupled with massless dynamical flavors. We will discuss its 
behavior, in a limiting regime of its parameters, at finite temperature and quark chemical 



potential (with both scales taken to be much below the UV Landau pole), by means of a dual 
charged black hole background (the case of extremely massive flavors has been analyzed in 
|2j). The gravity solution presented in this paper represents an improvement with respect 
to the ones found in the past by the authors and collaborators, in that it allows to scan an 
enlarged region of the (T, fj) plane. To be specific, our solution will be at first order in the 
flavor backreaction parameter e ~ \Nf/N c , but exact in the chemical potential. So, we will 
be formally allowed to explore the large charge regime of the system. 

As a first application, we will verify that the backreacted solution reproduces the thermo- 
dynamics calculated in the probe approximation [3] (in particular, the system turns out to be 
thermodynamically stable). This is a consistency check that the solution correctly describes 
the same physics as the probe whenever observables are computable in both frameworks. 
That is, the backreacted solution faithfully contains and extends the probe physics. 

Moreover, we will be able to check the consistency of the large charge regime in the above 
model as well as in more general "D3-D7 systems" (see the technical overview below). As 
we will see, in the parametrically large charge regime (including the zero temperature case) 1H 
higher corrections in the backreaction of the flavors, that is in e, are not subleading, spoiling 
the reliability of the solution at leading order in e. As a byproduct, this means that the probe 
approximation is not reliable in the IR at parametrically large charge or chemical potential. 
This fact is known, see e.g. [5], but possibly not fully appreciated. For example, it means 
that the use of probes to study the IR regime of the zero temperature, zero flavor mass limit 
of charged conformal systems is not trustworthy. Here we are able to confirm this fact from 
the explicit analysis of the brane backreaction. 

Finally, we will start investigating whether flavor dynamics affects the stability of the 
system. Our analysis extends previous results in the literature, where the flavors were treated 
as non-dynamical probes. As it was pointed out in several places (see e.g. [6]), there are 
various instabilities one could expect in our models. Charged fermions and scalars on the field 
theory side can give rise to various kinds of condensates: the corresponding symmetry broken 
phases should be possibly accounted for by novel dual gravity solutions. Instabilities could 
also show up in the spectrum of uncharged mesonic-like bound states. This is particularly 
interesting for us since it can be determined by studying (linear) fluctuations of dual fields 
around our charged black hole backgrounds. Focusing on a possibly critical scalar meson 
subsector, we will show that no instability appears at zero momentum in the regime where 
our gravity solution can be trusted. 

To help the reader going into the details of our analysis, we provide below a short technical 
overview on the holographic approach to flavored gauge theories, pointing towards the main 
issues raised in this paper. The reader familiar with the subject can safely skip (part of) this 
subsection. The way this paper is organized is then presented in a further subsection. 

1.1 Technical overview 

Holographically, in the ^ > 1, A > 1 limit, flavorless N = 4 SYM theory is described by 
type IIB supergravity in AdS§ x S 5 . The S 5 has a SO (6) isometry that relates directly to 
the symmetry of the scalar manifold in the field theory. Turning on temperature on the field 
theory side is accounted for by placing a black hole at the center of AdS*,. The lOd geometry 



'^With "parametrically large charge" we mean that the charge or chemical potential is larger by factors of 
N c than the other energy scales in the problems, as e.g. the temperature (this is the case explored in this 
paper), the flavor masses or possibly a dynamically generated IR scale (as in confining YM theories, such as 
the Sakai-Sugimoto model [3]). In particular, the extremely interesting case of zero temperature conformal 
charged systems with massless matter falls in this definition. 



originates as the near horizon solution of a stack of N c D3-branes. Introducing flavor into the 
theory is accounted for by adding a stack of suitably embedded Nj D7-branes [7]. The latter 
wrap an 5 3 C 5 5 , which translates into the 50(4) x U(1)r C 50(6) symmetry breaking 
pattern. The flavored SYM theory arises in turn as the low energy description of the open 
strings which end on the D3-D7 branes. Turning on a finite baryon charge density requires 
turning on an electric field on the D7 worldvolume. 

The D7-brane dynamics is described at leading order by the Dirac-Born-Infeld (DBI) 
and Wess-Zumino (WZ) actions. The finite value of the beta function, mentioned above, is 
accounted for by a running dilaton, which blows up at a radial position holographically dual 
to the field theory UV Landau pole. 

The treatment of the system consisting of type IIB supergravity with DBI and WZ source 
terms is generically a quite complicated problem to work out. We must, therefore, work in 
an appropriate approximation if we want to tackle the systemo 

Probe approximation. One of such approximations was started with the seminal paper of 
Karch and Katz [7], and goes under the name of probe (or quenched) approximation. It consists 
in considering a small number of fundamental degrees of freedom as compared to the adjoint 
ones: Nf/N c <C 1; then the 't Hooft coupling beta function is very small and can be effectively 
set to zero. This is the 't Hooft approximation. On the field theory side this corresponds to 
discarding effects of the dynamics of fundamental matter in the theory, which in the language 
of lattice gauge theory corresponds to quenching the quarks. On the geometrical side of 
the correspondence, the effect is to treat the D7-branes as probes on top of the undeformed 
AdS§ x 5 5 background. If these branes are all coincident, J\f = 2 supersymmetry is conserved 
and the flavor group is SU(Nf). Since in the undeformed background the dilaton is constant, 
there is no RG flow, signaling that in the field theory we have completely neglected the flavor 
field effects on the original vanishing beta function. As the gravitational content of the theory 
is kept frozen in the presence of the D7-branes, the dynamics of the flavors is described entirely 
by the DBI and WZ actions. 

Soon after the publication of [7], in which the super symmetric D7 probes were introduced, 
the spectrum of mesonic-like bound states in the supersymmetric theory was analyzed in [9]. 
In that paper, the spectrum of perturbations of several D7-brane fields, and their dual gauge 
theory operators, were studied and organized into the appropriate supersymmetric multiplets. 
Later, the theory was extended to include finite temperature physics [TUl [XT] and the presence 
of a finite quark density or chemical potential [3j [121 [13] . Following the construction of the 
thermal/charged backgrounds, perturbations of the system were studied in several regimes, 
see e.g. pH El EH E| • 

The contribution of the charged massless hypermultiplets to the free energy, to leading 
order in the probe approximation, was studied in [12], at finite or zero temperature. As a 
result the thermodynamics of the model was determined, showing that the system is thermo- 
dynamically stable. 

In |6j the possible occurrence of instabilities was studied by considering the mesonic spec- 
trum at finite charge and zero temperature for general flavor mass. In that paper it was 
outlined how the fluctuations of all the bosonic D7-brane fields effectively probe an AdS2 
near-horizon metric. This "effective" AdS2 region could be the remnant of a corresponding 
extremal region in the fully-backreacted background metric. In this case one would expect it 
to be related to a finite entropy degenerate system at T = as in the usual AdS-Reissner- 
Nordstrom (AdS-RN) solution. Such a degeneracy would then also be considered as a signal 



2 See [8] for examples of other solutions to this problem in the supersymmetric case. 



of possible instabilities (towards non-degenerate states) to occur. 

An especially interesting mesonic operator is the dimension-two scalar corresponding to 
a worldvolume scalar field of mass squared m 2 L 2 = —4, with L being the AdS§ radius. 
This mass sits precisely on top of the Breitenlohner-Freedman (BF) bound for AdS§, and 
a small perturbation could push the mode beyond unitarity. By studying the quasinormal 
modes of that scalar (corresponding to poles in the retarded two-point function of the dual 
mesonic operators) , the authors of |6j found that none of them corresponds to an instability. 
In turn, they found a "diffusive" mode (analogous to the "zero sound mode" appearing as 
a pole in the U(1)b charge density two-point function) with purely imaginary dispersion 
relation, Q = —iDq 2 , with q the momentum of the mode and diffusion constant D ~ 1/ ' \x 
from dimensional analysis. Near the origin of AdS$ spacetime, the mode is approximated by 
a massless scalar in AdS2- Since a zero mass scalar is above the AdS2 BF bound, the absence 
of instabilities consistently follows. In [18] it was shown that this mode is related to the Higgs 
branch of a moduli space of vacua (at zero temperature and finite chemical potential). 

In summary, despite the suggestions that the model should present an instability at low 
temperature, none could be found. Since all the relevant bosonic fields were included in the 
analysis of [6], it is possible that the instability can be seen only when backreaction of the 
D7-branes on the geometry is included. In view of the previous claim on the validity of the 
probe, however, one has to take in mind that the results described above are on solid grounds 
only away from the zero temperature, zero flavor mass limit. 

Smearing approximation. A different limiting regime, considered in the present paper, is 
provided by the Veneziano approximation, in which both Nf and N c are taken to be large with 
their ratio Nt/N c fixed. In this way, the dynamics of the quarks (in a perturbative language, 
the quark loops) is not ignored. In the holographic context this consists in considering the 
backreaction of the D7-branes in the original AdS§ x S 5 geometry. 

When the flavor branes are coincident, there are limitations (at least in flat space) in the 
maximum number we may have (see [19] for an explanation in this context). To avoid this and 
other (more technical than physical) problems, the D7-branes are distributed homogeneously 
in the two dimensions perpendicular to their worldvolume. This strategy was initiated in 
|20t I21| . The transverse distribution breaks part of the supersymmetry as well as the flavor 
group which reduces to its maximal Abelian subgroup U{l) N f . This is usually referred to as 
the smearing approximation. 

As already stated above, N = 4 SYM theory in the presence of fundamental matter de- 
velops a Landau pole at a certain energy scale (governed by Nf/N c ). Holographically this 
appears because the D7-brane backreaction sources the equation of motion for the dilaton, 
which does not admit a constant solution any longer. Actually, the dilaton shows a loga- 
rithmic divergence at a finite radial position r = r^p corresponding to the field theory UV 
Landau pole. In particular this implies that a conformal boundary is no longer available and 
the existence of a holographic dictionary is in dispute. However, experience with this and 
other systems suggests that physics of the IR can still be defined by applying the standard 
holographic dictionary at a large cutoff radial position below Tlp- This is the approach we 
follow in this paper, where we formally push the position r^p to infinity. 

D3-D7 systems in the smearing approximation. A N = 1 supersymmetric solution 
describing the backreacted intersection of D3 and smeared D7 branes at zero temperature 
and charge density was derived in [22| . The construction is quite general in that it is easily 
extended to infinite classes of flavored J\f = 1 quiver theories which arise when the D3-branes 
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Table 1: Five-dimensional fields originating from the ten-dimensional ones, taken from [28] 
G and A represent the lOd metric and world-volume gauge field, respectively. 



are placed at the tip of toric Calabi-Yau cones o 

The solutions in |22] are the starting point of the construction in |23j where their black 
hole (finite temperature) versions were found. Since the models are no longer supersymmetric, 
the equations of motion one needs to solve are second order and the analysis was restricted 
to small backreaction parameter e* ~ \*Nf/N c (where A* is the 't Hooft coupling at a given 
scale, see the text for details) to keep analyticity. The case of the charged black hole was 
later studied in |24[ [25] by including the effects of finite, small charge density per flavor 
perturbatively in the black hole solution of [23] (which is itself perturbative in e*)o In the 
present paper we generalize this work by giving the exact, in charge density per flavor, charged 
black hole solution (again, still perturbative in e*). As discussed before, we are able to use 
this solution to check the validity regime of the small e* approximation, and of the probe 
approximation as well. 

A reason to construct this solution is to study the effects of backreaction in the stabil- 
ity analysis of [6] recalled above. To perform this analysis we need to make a consistent 
perturbation of the IIB+DBI+WZ equations of motion. To do this we use extensively the 
results of [28], where a consistent reduction of the IIB+DBI+WZ action to five dimensions 
was presented. In this reduction the fields appearing are: 

- fields preserving the complex structure of the internal 5d manifold on which one com- 
pactifies (to ensure the consistency of the reduction) ; 

- fields neutral under the U(1)r R-symmetry, which is preserved by the smearing (this 
condition can be relaxed, see [291 EQl EI])- 

We reproduce table Q] from [28] where the field content of the 5d reduction and its lOd origin 
is specified. The action and equations of motion dictating the dynamics of these fields are 
written in [25] . We will not write them here explicitly but will refer to them extensively. The 
field strengths follow an equivalent nomenclature to the potentials shown in the table. For 
example, the lOd field strength three- form, F3, gives rise to three different 5d field strengths 
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i+ (a one- form), F% (a two- form) and F^ (a three- form). We refer again to [28] for exact 
definitions and details. 

The interesting scalar field sitting on the AdS§ BF bound (the one giving rise to the 
diffusive mode in [6j) corresponds in the former classification to the scalar .4.0 coming from 
the world-volume vector on the D7-branes. This scalar, in fact, is dual to an operator of 



In the flavorless case, when these theories are conformal, the dual gravity backgrounds have AdSs x Xse 
metric, where Xse is the 5d compact Sasaki-Einstein base of the given cone. 
4 These solutions have been studied further in [261 127] . 



dimension A = 2 + <5/|j where Qf is a measure of the backreaction (related to e* above by 
a factor of the dilaton, see (|2.8j) later). We could wonder whether the lim^o ^ = limit, 
found for the diffusive mode in the probe approximation, is perturbed in the backreacted 
case. Even if e* is taken to be small, it could weight positive or negative corrections to this 
limit, leading in one case or another to an instability (depending on the conventions: for us 
purely imaginary modes with Im(S7) > will be unstable, i.e. exponentially growing). 

1.2 Organization and main results of the paper 

This paper is organized as follows. In section [2] we present the minimal effective 5d gravity 
action relevant for the holographic study of the charged D3-D7 system with massless flavors. 
We will setup our perturbative approach presenting it in a slightly different (though physically 
equivalent) way as that followed by the authors and collaborators in the past (see ref. [28]). 
As we will see, this approach is quite useful in that it automatically identifies the fields 
contributing to the effective action at each order in the e* expansion. In particular it allows 
one to holographically study how irrelevant operators are "integrated" out along the RG-flow 
towards the IR. Moreover, as we will see, it shows how the effective action, order by order, 
assumes quite simple forms (e.g. Einstein-dilaton, Einstein-DBI) which are commonly used 
in bottom-up holographic setups. Finally, the approach allows to easily recover the already 
known uncharged perturbative black hole solution. The novel content of section [2] is a charged 
black hole solution at first order in e*. The difference with the solution presented in ref. [23J 
is that the parameter related to the charge density per flavor is treated in full generality in the 
present work, whereas in previous solutions it was treated perturbatively. This complicates 
the explicit form of the solution but increases its regime of validity. 

In section we study the thermodynamics of the charged solution at first order in e* and 
compare it to the results in [23J and with those obtained in the probe approximation. We will 
show, in turn, that the system is thermodynamically stable. 

In section [31 by studying second order corrections to the charged solution of section [21 we 
will show how in the perturbative-in-e* regime (and thus in particular in the probe approxi- 
mation with massless flavors) it is not sensible to approach the extremal T/fj, —> limit, for 
which an all-order solution is necessary. 

In section [5] we will study linear perturbations on top of the charged solution, focusing in 
particular on the mode dual to the operator with dimension A = 2 + Qf. We will find no 
instability in the regime of parameters where the solution is reliable. 

We present conclusions and final comments in section [6l Further useful results are collected 
in an appendix. 

2 The charged D3-D7 system 

A solution of the reduced system described in the introduction corresponds to D7-branes with 
a finite electric field in the radial direction. This describes holographically the presence of 
charge density per quark in the field theory. The study of this setup was initiated in [23] . 
where it is shown that the minimal set of fields that must be considered includes three scalars, 
<I>, w and /, dual to operators of dimension A$ = 4, A„ = 6 and Aj = 8 respectively, two 

(2) (2) 

vectors A\ and C\ and one two-form Cj , that upon redefinitions give rise to a massless 
vector field, a massive vector field, and a massive two-form, corresponding to a Aj = 3 flavor 



D For smeared anti-D7-branes this mode is not present, and Ao describes a A = 6 + Qf operator. 



current operator, and Ay = 2 + y/9 + cr 2 Qf and At = 6 + Qf operators]^ The parameter a is 
equal to — 1 (resp. to +1) if D7-brane (resp. anti-D7-branes) are introduced in the setup (see 
[28]). The case we will focus on is the a = — 1 one, though some of our results are presented 
for the general case. Schematically, the operators dual to the fields present in our solution 
are of the form [32l [331 1221 ES] 

0$ = trF 2 , O f = tiF 4 , O w e tr(W a W a ) , 

Of = tAN^ + iq ]m D tx q m - iD^q^ m q m , 

O v ,O t G TriWaWpW 13 ) + ... , 

where W is the gluino superfield, ip a doublet of spinors and q a doublet of squarks. 

2.1 The 5d effective action 

A 5d effective gravity action describing the dynamics of the minimal set of fields considered 
above can be obtained from a more general one |28| arising from a consistent Kaluza-Klein 
truncation of lOd supergravity. With the specific choice of the ansatz we will take below, the 
dynamics of the minimal set of fields relevant for the charged D3-D7 solutions is described by 
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l — f d 5 x v 7 - det g [R[g] +C s + C f + £ DBI ] , (2.1) 



C s = -f(df) 2 -20(dw) 2 -±(d$>) 2 -V($, f,w), 



V($,f,w) = 4eff +2w {e 10w -6)+-Q 2 eff- 8w+2 " + ^ef 
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C, = -^y- 8w (dC l + Q f dA l ) 2 -4e* + V +iw C 2 - I e *"f /(dC 2 ) 2 , 

Cdbi = -4Q /C * + ^ +!to ^l + i e -*-f/(tU 1 )2 . 

We have redefined 

C® = C 2 , C{ 2) = d , (2.2) 

and reabsorbed a 2na' factor in the definition of A.\. The constant Qf is related to the number 
Nf of D7-branes as follows 

«'■&■ (2 - 3) 

Here V(Xse) and V(X^) are the volumes of the five-dimensional internal Sasaki-Einstein 
manifold and of the three-cycle wrapped by the D7-branes, respectively^ 

The action above admits an AdS§ solution when Qf = and all the matter fields are 
trivial. The AdS radius is related to the constant Q c by L 4 = Q c /4. When choosing to work 
in L = 1 units (as we will mostly do here) one has then to set Q c = 4 accordingly|£| 

The 5d Newton constant G5 is obtained from the lOd one after reduction on the 5d internal 
compact space (of radius L) and it is thus given by 

1 L 5 V(X SE ) 



16vrG 5 " (2ir) 7 g 2 a 



a 



(2.4) 



6 A fourth scalar Cq corresponds to a Stiickelberg scalar coupled to the C{ vector and can be gauged 
away. The redefinitions leading to these dimensionalities can be found in [28] . 

7 In the flavored U = 4 SYM case, for example, V{X SE ) = V(S 5 ) = tt 3 ', V(X 3 ) = V(S a ) = 2ty 2 . 

8 Reinserting the correct dimensionalities one finds Q c = (g s a' N c (2tt) 4 )/V(Xse) where iV c is the number 
of D3-branes. 



Using the action ()2.ip . care has to be given to the fact that it does not strictly arise as a 
consistent truncation from lOd, but just as a "partially on-shell" reduction of the latter on the 
homogeneous "electric ansatz" we will choose. For example there will be extra conditions to 
be imposed on dC^ coming from the 5d equations of motion for H3 (which we take to be zero 
here). We will write these extra conditions explicitly, deriving them from the action given in 
[28], as an expression for the two- form field, i.e., we will not impose an ansatz for C2, since it 
will be determined from the other fields in our setup. 

The black hole metric ansatz we will consider is of the form 

ds\ = e 2A ^ [-b(r)dt 2 + dx t dxi] + e 2B ^ ^- . (2.5) 

Moreover, we will assume the scalar fields /, w, $ to be functions of the radial coordinate only. 
For the vector fields we will choose the following electric ansatz 

A 1 =A t {r)dt, d=C t {r)dt. (2.6) 

2.2 The perturbative approach 

Let us split the dilaton into a constant and a piece vanishing at the particular position r* 

$(r) = $* + <j>{r) , 0(r*) = O. (2.7) 

The thermal smeared D3-D7 systems studied so far (and in this paper) are perturbative in 
the parameter [23] 

e * = e Qf ~ i^v(x SE ) K ir c ' (2 - 8) 

where the 't Hooft coupling A* = 4mg a e *N C is defined at a certain scale set by radial co- 
ordinate r = r*. We can analogously define a running parameter e(r) = Qte®( r >. Since 
the flavor perturbation induced by the D7-branes and dual to the dilaton field is marginally 
irrelevant, we know (see [23]) that the shifted dilaton is logarithmically running to leading 
order: <p(r) ~ e* log (r/r*). This implies that, perturbatively, the beta function for e(r) is 
proportional to e 2 . This in turns implies that, at first order, the differences between e* and 
any other allowed value e(ro) are subleading in e*. When working with the first order thermal 
solutions, whose horizon radius r^ is holographically related to the field theory temperature, 
we will then be free to replace e* with eh = e(?7i) at first order in the backreaction parameter. 

As it was discussed in |23j . working in the e* <C 1 limit, allows us to decouple the IR physics 
(to which we are ultimately interested in) from the troubling UV Landau pole. Actually, the 
limit allows us to trust our IR description up to an arbitrary "large" radial cutoff r = r s 
(where we impose our solutions to match with the uncharged T = ones found in |22j). 
which has to be smaller than the one corresponding to the UV Landau pole scaleO 

Our approach to the perturbative solution here will be slightly different (though physically 
equivalent) to the one usually followed so far in the literature. The standard treatment consists 
in first writing down the exact-in-backreaction equations of motion and then expanding them 
in series of e*. In this paper, instead, we will first expand the action and only then derive the 
equations of motion. This approach is particularly useful if one is interested in determining 
holographically the IR properties of the D3-D7 plasma, like the thermodynamics. Since the 
field theory partition function is obtained from the on-shell (renormalized, Euclidean) 5d 



9 Actually, it must be r s < r a < tlp, where r a is the radial position where the holographic a— function 
presents a singularity |34| . 
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action, expanding the latter in powers of e* allows us to understand which fields will actually 
contribute to each order. 

Moreover, unless specified otherwise, we will formally send the cutoff r s to infinity, which 
physically amounts in neglecting all power-like r/r s terms in our solutions. This will allow us 
to automatically focus on the deep IR solutions which are those relevant for describing the 
thermodynamics. 

Before presenting the charged perturbative solution, let us first review what happens in 
the uncharged case using the approach described above. 

2.3 The uncharged case 

Setting the AdS radius to one, the effective 5d gravity action in the uncharged case reduces 
to 



S5 = iio;I d5x ^^~ 9 



R[9) - f(df) 2 - 20(dw) 2 - \{d>5>) 2 - V(*,f,w) 



(2.9) 



where 



V($,f,w)=4e 



>-f+2w / lOto 



6 + Q f e*) + lQ}eff- s ™+ 2 * + : 



>-f 



(2.10) 



When Qf = the action admits an AdS$ (black hole) solution of unit radius where $ =const 
and / = w = 0. The perturbative solution in e* = Q/e** <C 1, will arise in form of an 
expansion around the unflavored AdS§ background. In particular, we will expand the scalars 
as follows: 



<J>(r) = 


= $* 


+ e*Mr) + 0(el), 


f(r) = 


= e* 


-iW^ /i(r) 


w(r) = 


= e* 


f 1 Pf / X 

h \ — w-\ (r) 

60 V 40 V ; 



+ 0(ei 



+ 0(el 



(2.11) 



The constants in /, w are chosen in order to diagonalize the action and to have canonically 
normalized kinetic terms. Plugging these expressions in (|2.9|) and requiring that e*ipi <C 1 for 
ipi = (fi,wi,4>i), we find, to second order in e*, 



with 



and 



/ tfi t i 1 rlrt n 


R[g] + 12 (l - 


J~72) +€ * m 


167rG 5j / d V dCt9 


Cm = ~\{dh) 2 ~ \{dl 


oi?-\{d^) 2 


-v{<f>x, fi,wi), 



v(<p 1 ,f 1 ,w 1 ) = 4cj )1 + lGf 2 + Gwi. 



(2.12) 



(2.13) 



(2.14) 

This shows that the masses of the modes fi,Wi around AdS are m 2 r = 32 and m^ = 12 and 
thus (as already pointed out in the literature) that they are related to irrelevant operators of 
dimension Aj = 8 and A w = 6 respectively. The dilaton, instead, has no mass term and as 
such is dual to a marginal (actually a marginally irrelevant) operator. 
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The action to first order 



From the above expressions we see that, to first order in e*, the 5d effective action reduces to 
just an Einstein-Hilbert action with negative cosmological constant. This has thus a simple 
AdS§ (black hole) solution 



dsz 



r 
1? 



1 



3 



[-b (r)dt 2 + dxidxi] + L 2 ( 1 + 



dr 1 



3 7 r 2 bo(r) 



bo(r) 



1 



(2.15) 



with effective radius Lf given by L\ = L 2 [l + e*/3], where L is the ^dS 1 radius of the 
unflavored model. The thermodynamics to first order in e* can be thus easily obtained. All 
one needs is simply to replace L with Lf in the AdS black hole thermodynamical formulas. 



Thus, using T 



rh/(TrL 2 r) for the temperature, the entropy density s 



[l/(4G 5 )]7r 3 T 3 L 3 f as 



well as the other thermodynamical observables follow. The results are in perfect agreement 
with those found in |23j using the lOd solution 10 ! and with the results obtained in the probe 
approximation |11| . 



The action to second order 

At second order in e* the scalar fields start playing a role in the effective action. The equations 
of motion for f\ and w\ (dual to irrelevant operators) admit simple trivial solutions f\ = 
wi = 0. These precisely match with the solutions found in [23] when all the power-like cutoff- 
suppressed terms are neglected. In this way these fields are effectively integrated out and the 
constant values entering in the redefinitions (|2.1ip of /, w can thus be seen as the IR remnant 
of integrating out the irrelevant fields. 

The remaining non-trivial part of the effective 5d action, which only contains the metric 
and the scalar field <pi, can be rewritten, to second order in e„ asa simple Einstein-dilaton 
action of the Chamblin-Reall 1 35 1 kind 



1 f d 5 r 


y '- det g 


m - \(d^) 2 - 


- F e w 


16vrG 5 J d X 


the order we are interested in, 




Vq = -12 


1-^- 
3 


72 _ 


e* 

' 7 = y 





S *ff = TTTV^ / d'xy^detg R[g] - -(dip) 2 - V e^ , (2.16) 



(2.17) 



Despite the fact that we have to stop at second order in e*, let us recall, reporting e.g. the 
results in [36], that this action has a known exact black hole solution. In the r = <p "gauge" 

it reads 

dp 2 



dsc 



3 2A( V )r 



-h(ip)dt + dxidx^ + e 



2B{y)_ 



h{fp) 



(2.18) 



where 



„2A 



exp 



2 ip 

"3^ 



„2B 



37 2 , , 

2 — exp(-79?) 



67 2 ^o 



1 — exp 



8-3 7 2 
67 



(<Ph ~ <P) 



(2.19) 



10 Notice that the standard AdS radial coordinate used in (|2.15p and the r-coordinate used in [23] are simply 
related by a rescaling at first order: r hcrc = [1 + (5e„/24)]r t herc- 
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Prom these formulae we get the entropy density and temperature of the black hole as (see 
also [37]) 



,3Ai, 



fh 



4G 5 4G 5 ' 

pAu—Bh 

\h'{<Ph)\ 



4tt 



1 

47T 



:Vn\ 1 



3 ,2^-|^(2-37 2 ) 



t7 e 



so that, to second order, 

where 
so that 



8 = a T 3 



1 4- ^ 4- — f 2 
1+ 2 + 2A Eh 



eh = e* 4 el log(7rT) 



T deft - J 



Above, <jqT 3 is the entropy density of the unfiavored plasma, with 



3r3 



n 6 L 



tt 5 N? 



O-Q 



4G 5 2F(X 5S ) 



(2.20) 

(2.21) 

(2.22) 

(2.23) 

(2.24) 



being a measure of the number of degrees of freedom of the unfiavored theory (it is proportional 
to the central charge holographically given by a = N 2 7r s /(4y(XsE))- The entropy density, as 
well as the "RG-running" formula for e^ given above, precisely coincide with the ones found 
using the full lOd action in [23]. The remaining thermodynamical observables also follow 
accordingly. 

These results confirm that in the deep IR limit (i.e. for T <C A s <C Ajjv) the thermody- 
namics is captured by just the marginally irrelevant operator dual to the dilaton field. The 
same conclusion holds for the hydrodynamic behavior, as it has been shown in [26] . 



2.4 The charged black hole solution 

Let us consider now the charged case, which is the focus of the present paper. In the following 
we will adopt the same redefinitions for the scalar fields as in (|2.1ip together with the following 
ones for the forms 



Ai = e^' 2 A , dC 2 = Q f e**/ 2 F 3 , d = Q f e^' 2 V . 
With these redefinitions we find, at second order in e* 



'e// 



16vrG, 



— / dxV-dets 



%] + 12- Ir :V /l + ^+ e 2£ 2 



(2.25) 



(2.26) 



where F = dA, Y = dV and Fa = dC 2 - The term C 2 is given by 

(dhf (d Wl f (dfa) 2 5 , 2 



£2 



2 2 2 6 3 L J 

-Gi[F 2 ]/i " tfil^Vi - I6/1 2 - 6w 2 - 40xA![F 2 ] + 

-\(Y + F) 2 - l -F 2 -W 2 , 



(2.27) 
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with the functionals 



Ai[F 2 ] = 


x/1 + C^/2) ' 


d[F 2 ] = 
HAF 2 } = 


3V2F^ + 16V2 - 16V2 + F 2 


V15V2 + F 2 
-2\/2 + 2V2 + F 2 



(2.28) 
VD 

In the uncharged case F = V = F3 = we get Ai = 1, G\ = H\ = so that the uncharged 
second order effective action (|2.12p is recovered. 

The action to first order 

From eq. (|2.26p we see that, at first order in e*, the effective action reduces to a simple 
Einstein-DBI one, where only the metric and the U(1)b field F enter. The thermodynamics 
of the charged D3-D7 system at first order in e* is thus fully captured by that simple action. 
Using the electric ansatz given in (|2.5p . (|2,6p we find that the equation of motion for the 
vector field is readily solved by 

A(r)+B(r) r 3 

F r t = , d , (2.29) 

/ e 6A(r) + r 6 



where r^ is a dimensionful constant naturally associated with the charge density per flavor rH 
Writing the metric as ds 2 = gMN( r )dx M dx N it is easy to see that Einstein's equations 12 ! 

2 
Rtt + ^9tt + ^ e " 



9tt J r 6 + „3 
3/2 V d + 9xx 
Qxx 


-2e*- 

4 « 

-2e*- 


9tt9xx 


= n, 


\J r d + 9 XX 




*<-xx ~r ^9xx 


\J r d + 9 XX 


= 0, 


1/2 
9xx 

3/2 
9rr9xx 


9rr . / r 6 4. n 3 
3/2 V d xx 
9xx 


= n 


J r d + 5W 





*rr + *9rr + fe.^ \l ^ + ;#,. 2r, ;^TI_ . (2.30) 

admit the following simple solutions (found imposing regularity at the horizon and UV match- 
ing with the uncharged flavored solution in (|2.15p ) 



9* = ~r 2 (l - |) 6(r) , 
Ixx = r 2 (l - - 
.(! + ¥)' 



</ 



r 2 b{r) 



"Notice that rd is not a radial position, since it can have a negative value in correspondence to the sign of 
the charge density per quark. 
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Cfr. e.g. [35], eqns. (18)-(20), setting V = 0, Z x = Z 2 = 1, A = 1, <j> = 0, d = 4, A = -6,T b = 4e* 



14 



where 



b(r) 

b q (r) 



1 + 

1 
3 r 3 



c* 



A 



+ e*6 q (r) 



rD + r ° + ^W r ' + r ° 



rjG F (r)-G F (r h ) 



1 / r 6 11 

3 \r + ^ 6 3 



3V4 



F 



cos 



(1 - yg) r 2 + r 2 
(1 + x/3) r 2 + r 2 



y/3 



(2.31) 



-F[x|y] is the Elliptic integral of the first kind, B(x; a, b) the incomplete beta function and we 
have chosen the integration constants so that &(r/J = 0. Finally, the electric field reada 13 ! 



F 



rt 



r D +r b 



(2.32) 



This result is in agreement with the solution obtained in the probe approximation in [12] (see 
their eq. (3.8)). Integrating the expression above and imposing At{rh) = we find 



Mr) 



rd 



G F (r) - G F (r h ) 



(2.33) 



The action to second order 

Let us now consider the effective action to second order. If one is interested in the thermody- 
namics (obtained from the free energy, which in turn is holographically related to the on-shell 
gravity action) , it is clear that the relevant equations of motion for f\ , w\ , C2, (pi and V which 
appear in the action at 0(e 2 ), can be solved just by computing them on the "zeroth-order" 
background 



dsf 



-bo(r)dt + dxidxi] + 



dr 2 



r 2 bo(r) 



b (r) 



1-^ 



& 



rt)0 



ty*\) I fy*C 



, so that (-F 2 )o 



r o _|_ r > 



6 ' 



(2.34) 



The equations for /1, w\, (pi, V and Ci which we need to solve are 



(5r 4 -^)/{+r(r 4 -^)/{' 

(5r — r h )w'i + r(r — r h )w'{ 

(5r 4 -r£)# + r(r 4 -rM 

d-k F 3 

d* (Y + F ) 



r 3 (G 1 [(F 2 ) ]+32f 1 ), 

r 3 (i/ 1 [(F 2 ) ] + 12u;i) 

4r 3 Ai[F 2 ], 

0, 

8*0^, 



(2.35) 



where the Hodge dual is computed on the zero-th order background. It F = these equations 
consistently reduce to the uncharged ones, which have /1 = W\ = as solutions. In the charged 
case these trivial solutions are not admitted and all the scalar fields are running. As in the 
last equation F = Fq, we see that the field V (a massive vector field dual to an irrelevant 
operator) is "effectively" decoupled from the other fluctuations in this setting. 



In order to get the first order correction to this expression we need to consider the second order effective 
action. See below. 
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The above system of equations is actually not complete: there is in fact a further constraint 
which should be imposed by hand which comes from the requirement that the 5d ^3 fields 
(which do not appear in our effective 5d action) can be consistently set to zero. In particular, 
from eq. (78) in [28] we get that to leading order 



F 3 {r) 



r 3 V ^ F ^ H^dx 1 A dx 2 A dx 3 = r 3 d dx 1 A dx 2 A dx 3 , 



(2.36) 



x/1 + (^ 2 )o/2 
which readily solves the corresponding equation in (|2.35p . 

The solutions to the rest of the equations (|2.35p (found by requiring regularity at the 
horizon and UV matching 14 ! with the T = uncharged flavored solution of section I2.3[) are 
quite involved and we can give them just in a semi-analytic form 



-Jl 



r 



h 



1 1 2r 4 
"40 ~ 902rJ-r' h 

+ r 4 r 4 (2r 4 -r 4 ) 



(2r 4 



r 4 

r h 



f 

J r 



/V«0 _l_ /V»U 



40(2f 4 - r\) 2 
G F {f) - G F (r h ) 



d:f 



(2.37) 



20f(f 4 -r 4 )(2f 4 -r 4 ) 2 



dr 



-W\ 



1 , p v 2 [ 2 k 

60 M 2 ^f 

r6R ^ 2 ( 2 I - 

10r 4 



1 



1 
60 



Pi/2[2%~1) r ?»R 1/2 [2^-\ 



10r 4 



f b + r° 



dr 



(2.38) 



Ri ' 2 \ 2 k ~ 

Pl/2 (2^J 



r b P 1/2 ( 2 



7Z 



rh 



f 6 + r*\ 



dr 



P1/2 2 



^ 



1 



r s f 6 P 1/2 ( 2 



^ 



r" + r„ 



dr 



where we have used the explicitly real combination Ri/ 2 (x) = Qi/2( x ) + ^; Pi/2i x ) ■> with P n (x) 
and Q n (x) the Legendre functions of the first and second kind. 
The dilaton solution reads 



T\ r° + r~. 



r h \/rl + r% 



K 



log b (r) 



r 6 + r 6 



r h\/ri + r 



K 



r 4 1 l' r * 2f 6 -I- 

+ ^{D F {r)-D F {n)) + —, I 



: log 6q (r)dr . 




where we have defined 

D F (r) = logb (r)(G F (r)-G F (r h )) , 

and finally, for the vector field, we write 

4 4 



(2.40) 
(2.41) 



so that the equation of motion can be written aa 15 l 



8 r 



±J{ 



8r- 3 J x 



u 



4 / r 6 +r t 



14 We present the solution at finite UV cutoff r s . 

15 Recall that the parameter a is equal to —1 in our setup and to 1 in the anti-D7 brane case. 



(2.42) 
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with solution 



Ji=a{2 



r 2 + ( 1 + r 4 

'h V ' h 



log 



r Z _|_ r i 



((V3-l)rg + 2rg 

256r2 r 2 



-(G F (r)-G F (r h )) (2-43) 



+ 



;i+>/3) 



r b + r c 



128rM ( (l + V3)r2 + r 2) (,* + (l + >/5) rg) 



128r2 



(G s (r) - G E {r h )) 



a 



4 i 4 
r + r ft 



2L i + I ! + 4 I lo ? 



((V3- lM + 2rf) 



256r2 r 2 



(G F (r s ) - G F (r h )) 



+ 



! / (1 + yg) r^rg + j (l + yg r^/j + r£ 
+ (1 + VS) r?) (rj + (1 + V3) ri) 

(f 4 + r 4) 



128r* 

V3r| 
128r? 



(G E (r s ) - G £ (r h )) 



+ a- 



r + r^ 
Ur & h 



log 



tv*\j I /V" *-■ 



(7 

96 



1 + 



log 



r d + \ r b + rj 



rj + \ r b , + r\ 



r l + r s 



r l + ri 



df 



where Ge is defined in the same way as Gf in (|2.3ip . but with the elliptic function of 
the second kind i£[a;|y] instead of the first kind one .F[x|y]. In particular, notice the linear 
dependence with the WZ factor a, describing the dependence of this function on the charge 
of the smeared branes. 

While for the other fields the zero-th order background suffices, in order to solve for F r % we 
need the first order background metric (|2.3ip . which has gttdrr = — 1 and y/g = r 3 [l — (e*/2)]. 
Dubbing R d the integration constant, we find that, to first order 



Frt(r) 



R 3 d 



l + e* 



r e + R 6 d 



QM 



(2.44) 



, r 6 + i? 6 

where Q[r] is expressed in terms of the on-shell values of the scalar fields and the vector V 



QM 



r 3 <9£ 2 
~2RjdF7t l °. 



(2.45) 



and the derivative has to be evaluated on the zero-th order solution for F r t. 
The integration constant R^ can be written as an expansion 



Rd = r d [l + —pi 
with p\ <C 3/e*. This way we can equivalently write 



(2.46) 



Frtir) 



r b + r o 



1 + e*- 



tyKi I /y' 



:(Pi+Q[r]) 



(2.47) 
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with (to be explicit) 



Q[r] = g 



1 + 



-l + 60i)r§ , h(2r 6 + 5r 6 d ) 



6r 6 



15r 6 



-Tl* 




(2.48) 



Comments 

The whole charged solution described in this section is perturbative in e* but exact in r<f. In 
|24j the solution with |r*d| <C r/i was presented, which implied a second perturbative expansion 
in S = r\/r\. The solution in [23] is thus only valid when the energy scale of the charge density 
per flavor (or chemical potential if a change of ensemble is performed) is much lower than 
the energy scale associated to the temperature. The fact that our solution is exact in r^ does 
not mean that the solution is valid for asymptotically large values of the charge density (or 
chemical potential) per flavor, since the solution is still perturbative in e* and large values of 
rd may push the solution out of the region of validity. This will be clarified in the following 
sections by studying the extremal limit. This limit requires the knowledge of the solution 
at second order in e*. This solution is more involved and we have not been able to find 
a complete analytic expression for the metric in that case (the explicit solutions to f\, w\, 
4>i, etc. enter in the equations of motion, implying that an analytic approach becomes an 
extremely convoluted task), however we can obtain crucial informations that allow us to find 
some conclusive results. 



3 Thermodynamics 

In this section we discuss the thermodynamics of the charged solution at first order in e* ~ e/j , 
where we have an analytic solution for the 5d metric and the electric U{\)b fieldl 16 ! 

To begin with, let us notice that the thermodynamics of the same system has been cal- 
culated in the probe approximation in [3j. Let us immediately make it clear that, in the 
cases we can compare, the thermodynamical quantities calculated here and in |3j will coincide 
modulo the temperature. That is, the observables will coincide once we write the result in [3], 
where the temperature was the constant temperature of the unflavored theory, in terms of 
the temperature of the backreacted solution] 17 ! 

The temperature of the charged black hole solution found in the previous section reads 



T = — 

IT 




T n 




(3.1) 



where To is the temperature in the unflavored case. Equation (|3.ip is a monotonic function 
of r^ (namely, of the chemical potential or charge density per flavor) vanishing at 



extr 
r d 



4 



1/6 



Vh- 



(3.2) 



We do not expect, however, that we can really reach the extremality regime where T ~ since 
there the corrections to Tq are order unity and we are not guaranteed that higher orders in e^ 



^Extrapolations to QGP RHIC temperatures give rise to en ~ 0.24 |23j and we will use this value where 
needed in our plots. 

17 This was the case also in the perturbative solution in |24j . 
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will be subleading. In fact, we will argue in section [4] that the corrections are not subleading, 
so that the extremal perturbative solution is not reliable. 

The field theory entropy density, s, holographically identified with the Bekenstein-Hawking 
entropy of the dual black hole 

(3.3) 



can be recast in terms of the temperature as 



a T 6 



where 



4g/* I 1 " 


2) ' 


re as 




l-^ + e.Vl + 5 2 


6=^ 





(3.4) 



(3.5) 



and co defined as in (|2.24p . It is easy to see that the entropy density precisely reduces to that 
found in [21] in the (5<1 limit at second order in 5. 

The charge density n q is given by the holographic relation 



11 a 



L 5 dA t UttG 5 L 5 



so that 



■;?.,, 



IT' 



T3~ € h\ 






V(X 3 



7T 



7/2 



16vr V(X SE ) 3/2 



X h N f N c 5. 



(3.6) 



(3.7) 



This relation matches (to first order in e^) with that obtained in [24] in the small 5 limit. This 
is related to the fact that (no matter the value of 5) at this order the field At is effectively 
decoupled from the vector and two- form fields in the action (|2.1|) . Notice that the baryon 
charge density ub is related to n q by the simple relation n q = N c hb- 

From (|3.7|) we get that at first order in e^ and fixed charge density (canonical ensemble) 






36 



(3.8) 



Using this result it is easy to verify that the Helmholtz free energy density, given in terms of 
the incomplete beta function 



£(<*) 



1 



4 

£4/3 



croT 4 



B\\ 



_ £h 

2 

1 1 

6'3 



1 - 2 Vl + £ 2 + £(5) 



B 



1 



1 1 



l + 5 2 '6'3 



(3.9) 



satisfies the thermodynamical relation s = —df/dT (canonical ensemble). Notice that £(5) ~ 
3<5 2 in the 6 — > limit. Using this result we see that the formula above precisely reduces to 
the corresponding one found in [21] in the small 5 limit. 

The energy density (holographically related to the ADM black hole one) reads 



a T 4 



1 



ei, 



1. 



i_ 2 yr+^-|£(<5) 



(3.10) 



and can be deduced from the relation e = / + sT. From this expression we easily determine 
the pressure p, and thus the Gibbs free energy density u, as 



P 



-ui 



+ 0(e 2 h ) 



(3.11) 
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0.5 1.5 2.5 

Figure 1: 3/3^=4 (solid line) and e/e^=4 (dashed line) as functions oird/rn with e^ = 0.24. 



This also follows from the fact that at first order in e^ the trace of the stress energy tensor at 
thermodynamical equilibrium does not receive corrections with respect to its (zero) value in 
the conformal limit. In fact we know that the conformality breaking effects are higher orders 
in e/j |23j and that the chemical potential does not contribute to the conformality breaking 
in a theory like the present one where matter fields are massless. This also implies that the 
speed of sound is v 2 s = 1/3 and the specific heat at fixed chemical potential is cu = 3s. 

In figure [T] we focus on the flavored N = 4 SYM case and we report the entropy density 
and the energy density, normalized with respect to the ones in the unflavored theory. Nothing 
special happens at the would-be "extremality point" r^ xtr ~ 2.32r^. In particular, the entropy 
density is finite. 

As said above, in order to compare with the probe results in [S\, we have to normalize our 
quantities by the appropriate power of T/Tj^ = 4. After this normalization, the agreement of 
our results with the ones in [3] is perfect. 

We can check that the thermodynamical relation s = —duj/dT at fixed chemical potential 
(gran-canonical ensemble) is consistently satisfied . First of all notice that from the thermo- 
dynamical relation \xn q = f — to and from formulas (|3.6p . (|3.7p we get (to leading order) 



T 



It can be checked that, consistently, 



7T E(<?) 
12 5 



irX h 



V{X, 



SE) 



(3.12) 



=**/2 



/' 



27ra' 



7 A t (oo) , 



(3.13) 



where At{r) is given in (|2.33p . In the small 5 limit these expressions precisely agree (to first 
order in e/J with those given in [24] . In particular in this limit we have [ijT ~ 5 hence [ijT 
scales like r\. In the opposite S S> 1 limit (which can be formally taken with the proviso of 
taking r^ not greater than r^ xt ) we get 



7T 



ttA/j 



// 



T 12 V V(Xse) 



<4i> ,/3 - 6 



TT 



vrAft 



12 W V(X SE ) 






(5»1), (3.14) 
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Figure 2: ju/v / AT)v=4 as a function of r^/V^ in the perturbative charged (dashed line) and 
exact charged (solid line) cases, with eh = 0.24. 



which shows that fi/T scales linearly with r^. This behavior is made evident in figure [2] where 
(focusing on the M = 4 SYM case) we also compare our results against the perturbative-in-e) 
case, normalized with vA^TV=4- From (|3.12p we get that 



dTj 






(3.15) 



which is what we need to verify that s = —dco/dT at fixed /i. 

Now we are in the position to calculate the susceptibility matrix, in order to check the 
thermodynamical stability of the system. The "quark" susceptibility \ = d^n q and xtt = 
— 752^f are reported in figure [3] for the flavored J\f = 4 SYM master example. As we can see, 
both susceptibilities are increasing positive functions of the chemical potential. Obviously, 
their values at \x = are the same as in the perturbative case 



Xpert 



4vr z 



e h N 2 c T 



XTT,Af=4 



3vr 2 



N^T 2 



(3.16) 



The off-diagonal susceptibilities are as in the perturbative case 

.2 /_ \3 



XTp, = X m t = d T n q = -^jthN 2 [^-) T 2 . 



(3.17) 



Actually, we do not need their expression at leading order in e^ for the issue of stability 
(i.e. the positivity of the determinant of the susceptibility matrix). The point is that these 
observables are already of order e^, so they are subleading in the determinant, which at leading 
order results to be simply proportional to the "quark" susceptibility % (times the unflavored 
part of xtt, i-e. xtt,N=a)- 

In conclusion, the thermodynamical stability of the system is guaranteed by the positivity 
of the "quark" susceptibility % f° r every charge value. 
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Figure 3: x/Xpert (solid line) and Xtt/xtt,M=a (dashed line) as functions of r^/r^ with 
e h = 0.24. 



4 The issue of extremality 

As we will see in the following, the charged solution in section [2] is not complete around T ~ 0, 
forbidding us to trust it at extremality. The point is that going to T = would require, strictly 
speaking, going beyond the perturbative-in-e regime. In fact, at first order, we have found 
that T = Tq[1 — 0(e/i)] and implementing the extremality condition (|3.2j) would imply taking 
the subleading 0{th) terms to be of order one. The perturbative approach would have some 
hope to be sensible only if, despite being at extremality, the second order terms keep staying 
subleading with respect to the first order ones. As we will show in the following, this, not 
unexpectedly, is not the case. 

The extremality condition for the charged solution at first order can be rewritten as 



y 



(n 



3a/ 1 



9 



f + 0(4). 



(4.1) 



where y = S^ 1 = r?/r5. In this way we can say that turning on the flavor backreaction, the 
extremality condition y = turns into y = e^/ '3 at first order in e/j. 

It is easy to show that taking this "first order extremality limit" the thermodynamic 
relations found in the previous section precisely reduce to those found at T = in the 
charged, massless, D7 probe approximation [T2]. In particular the extremal metric develops 



an ^(iS^ near-horizon region and the entropy density s 



n^/ y^ 



is non-zero. The free energy 



scales like rj accordingly. 

Now, in view of the previous observations, the question is whether this first order extremal- 
ity limit makes sense. To answer this question we may ask how the extremality condition is 
modified at second order. For dimensionality reasons (there are no other scales in the deep 
IR apart from r^ and r^) the T = condition to second order should read as 



y 



eft 



^^TTy^-elC[y}=0. 



Now, if C[y] is an analytic function of y, we can expand it around zero getting 

y = ? + 4C[0] + 0(e 3 h ) . 



(4.2) 



(4.3) 
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In this case, the second order corrections would be subleading and the first order extremality 
limit would not be spoiled. 

If instead, say, C[y] ~ ky~ m with m > when y — )■ the situation would be very different. 
Take, for example, m = 1. In this case the first order value y = e/,/3 would be modified to 
y = (e/j/6)(l + y/1 + 36k). This would mean that the perturbative expansion in eh would 
not be reliable when going to T = 0: higher order corrections to the T = condition would 
modify the critical value of y to lower orders. For m = 2 we would get that the critical value 
of y gets a zero-th order correction. This, again, would tell us that the perturbative expansion 
in €h is not reliable when going to T = 0. 

In the following, examining some properties of our charged solution to second order, we will 
find that, unfortunately, the above situation (actually the m = 1 case) is precisely realized. 

For our purposes, since we are just interested in the T = condition, it suffices to consider 
the near- horizon behavior of the matter fields in the small y limit, with y = r^/r^. From the 
analytic expressions reported in section [2j we find that 

V t (r h ) = 0, Y rt (r h ) = V t \r h ) = JJ(tt + log4 - 1) + 0(y 2 ) ~ 2.12 + 0(y 2 ), (4.4) 

t < \ 4 + 37r-181og(2)l 0.06 

fl(r h ) = t= &v - + ■■■ + ... (4.5) 



1>/T5 y 



w 1 (r h ) = —= + ...-0.1 + ... (4.6) 



and 



From these we see that, among these modes, the leading contribution to the IR physics 
in the extremality regime comes from the field /i. The horizon contribution of wi, as well 
as that of the field strength Y of the massive vector mode V, is a constant. It would be 
interesting to precisely understand why a different role is played by these fields, despite the 
fact that they are all dual to irrelevant deformations. 

As for the dilaton, its value at the horizon can be reabsorbed by a redefinition of e. In 
fact, as we know, eh = Qfe®( Th > = e*(l + e„,0i(r/ l )). 

The horizon value of F r t from (|2.47|) turns out to be 



F r t(r h ) 



1 



V^+i 



y 2 
y z + 1 



(4.7) 



and in the y — > limit 



<3M = ^T73 + 0(l/y 2 )^- + 0(l/y 2 ), a ~ -0.04, (4i 



were we have used the fact that at first order e*^»i(r > ^) ~ 0. From this expression we see that 
the T = extremality condition y ~ e*/3 which emerges at first order cannot be consistent 
with the perturbative expansion to second order. In fact, if y ~ e*/3 the C(e*) correction in 
the expression for F 2 Jrh) turns out to be not parametrically suppressed: it is just a number 
equal to 3a ~ — 0.12^£| This should be sufficient to state that the second order correction to 
the T = condition cannot be made parametrically subleading with respect to the first order 
one. 

To strengthen the above statement let us proceed as follows. Defining the off-shell matter 
Lagrangian as 

L m = -4e* y/l + F 2 /2 + elC 2 , (4.9) 



'Recall that by construction pi <C 3/e* ~ l/y. 

23 



we know that Einstein's equation are 

Rmn — -z9MN — Qqmn = -QMNL m — MN . (4-10) 

In order to get, from these, informations on the T = condition, it suffices to focus on the tt 
component. Our metric ansatz will be of the generic form written in (|2.5p and, as usual, we 
will require that b{ru) = 0, so that b(r) = b'{rh){r — r^) + 0((r — r/J 2 ) close to the horizon. 
The T = condition is realized when 

b'(r h )=0. (4.11) 

Now, taking the tt component of Einstein's equation, multiplied overall by g , expanding 
around the horizon, taking just the leading zero-th order terms (in an expansion in r — r/J, 
and using the on-shell horizon values of the various matter fields, we get that the T = 
condition is, as expected, precisely of the form written in (|4.2p . 

Moreover, we also discover that, in the y — Y limit, the function C[y] goes as 

C[y] = (4 + 3^-181og2)(l + 2vr-121og2) + ___ _ k + __ 

I20y '" y 

so that we fit in the m = 1 case discussed above. 

These results confirm that in order to reach the T = regime in the charged D3-D7 
models with massless flavors we need to abandon the perturbative-in-e approach. This implies 
in particular that the extremality limit cannot be taken in the probe approximation. This is 
something which has already been stated in [5] but perhaps not fully appreciated in the 
literature. The reason for the breakdown of the reliability of the probe approximation is 
simple: if T = the energy density of the unflavored conformal theory (say, N = 4 SYM) is 
zero. We cannot thus work in the probe approximation consistently, since even a very small 
number of charged flavors would provide a larger contribution to the energy density. For the 
same reason, in the case of theories displaying a dynamical scale A/#, the reliability of the 
probe approximation could be obtained only for a limited range of values of fi/AjR. 

A final comment is in order. The DBI action we have been using constitutes an approxi- 
mation to the full dynamical flavor effects. In particular, it re-sums the "one-window graphs" 
in the Veneziano limit, that is the graphs with one quark loop. At order e 2 one generically 
expects further corrections, e.g. two-window graphs, to be present. Moreover, despite the 
fact that they are unknown in the D7-brane case, there could be "thermal" corrections (given 
by the blackfold approach [39]) taking into account the proper thermalization of the brane 
degrees of freedom. In order for these corrections to spoil the result of this section, they 
should all conspire to cancel the non-subleading e 2 contribution discussed above. While we 
cannot discard this possibility, we find it very unlikely (for example, we suspect that in the 
T <S Auv regime we work in, thermal corrections would be suppressed); were it realized, the 
solutions presented in section [2] would include bona-fide extremal charged black branes, with 
AdS2 x R 3 near horizon geometry. 

5 Perturbations 

We consider now linear fluctuations of the 5d fields around the charged solution. We refer 
again the reader to [2H] for a complete list of the equations of motion of the 5d system, from 
where the equations of motion of the fluctuations can be derived by linearizing around the 
setup described in section [2j 
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EOM[A ] = 


= 0, 


EOM[\{[ = 


= e* <7[A 


£OM[Ai] = 


= o-[Ai] , 



Representing generic fluctuations of type IIB supergravity fields by A and fluctuations 
of D7 worldvolume fields by A, we can expand in powers of the backreaction parameter as 

A = Ao + e* Ai + • • • and A = e* Ai H where in the last expansion we have started at order 

0(e*) since in the absence of D7-branes there are no worldvolume fields. In an e* expansion 
the different fluctuations equations of motion can be written in the hierarchic way 

(5.1) 
(5.2) 
(5.3) 



where EOM[-\ represents linear, possibly coupled, second order differential operators and &[•] 
some linear, at most first order, differential operators acting as a source! 19 ! 

If we had considered a D7 probe with a trivial embedding in our setup, the equations 
governing the dynamics of this extra D7 would be given by (|5.2p with a [Ao] = (because 
the supergravity fields remain frozen in the probe approximation). In the present setup 
supergravity and worldvolume fields fluctuate with amplitudes of the same order and their 
equations are naturally coupled. 

We will consider only leading order effects, since the higher order in e* terms will contribute 
just with small corrections that will not change our conclusions. Working at leading order 
means that equations (J5.3J) . etc. are not considered. This implies particularly that the 
backreaction corrections described in section [2] do not appear explicitly in the equations we 
solve. 

5.1 Restricting the fluctuating modes 

(2) 

Due to the presence of non-trivial components in the background RR potentials C\ and 

(2) 

C% i a complete analysis of all the perturbations becomes very involved. However, we can 
consider a consistent truncation of such perturbation modes, which includes in particular the 
field we are interested in, Aq. We will consider that any field \P in the background solution 20 ! 
is perturbed by 

* -»• $ + e ikx 8^ , (5.4) 

where we can use the residual SO(3) little group invariance of the solution we are perturbing 
to set k ■ x = —uit + qx 3 . Since we will consider fluctuations in a black hole background, 
the related frequencies associated to the 5^/ modes will generically be complex. We will also 
have to pick up appropriate boundary conditions for the fluctuating fields at the horizon: 
incoming-wave boundary conditions at the horizon are known to be relevant for computing 
retarded correlators in the dual theory . 

We will restrict our analysis to the q = case in which the little group remains unbroken. 
As a result, the fluctuating modes classify into tensorial, vectorial and scalar SO(3) modes. 
As we explained in the introduction, we are mostly interested in analyzing the behavior of 
the fluctuation of the Ao mode, which is a scalar under 50(3). Therefore, we will focus just 
on the scalar fluctuations of the fields. These include (see also Table [1]) perturbations of all 



19 Notice that EOM[-] and <r[-] are not necessarily the same in every equation; we have omitted extra labels 
for the sake of clarity. 

20 In particular ^ can have tensorial indices which we are not writing, and can have the background value 
# = 0. 
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the scalars present in the system, and components of the vector and tensor fields spanning 
the t — r subspace. Components of fluctuations of the metric in the t — r directions, as well 
as its trace, must be considered as well. Counting the number of modes we have 

1x8 (scalars) + 2x5 (vectors) + 1x2 (2- forms) + 4 (metric) = 24 modes. 

This is still a large number of modes to analyze. We can make a further reduction on the 
number of fields by imposing that the equations of motion that do not vanish identically in the 
background (the ones for the three scalars /, w and <3?, Einstein equations and the equations 

(3) (3) 

for H^ and -F 2 ) do not receive corrections at first order in fluctuations. This implies that 
we will not consider fluctuations of non-vanishing background scalar fields, nor the metric. 

(2) 

Furthermore, consistency requires that fluctuations of the two-form mode B^ to cancel. 
Doing the counting we observe that now we have four scalars (SB^ , 5Cq , 5Cq and 5Ao) 
and the t and r components of three vectors (5B[ , 5C\ and SA\), adding up to ten different 
modes which form a closed system (in principle we could have considered the tr component 
of the two- form 5C2 as well, but it is pure gauge in this setup). In appendix [A] we list 
the equations of motion describing these perturbations at leading order, along with technical 
comments. From now on we skip the 5 in the names of perturbation fields. 

A schematic description of the operators dual to the fluctuating fields we consider is: 

- the axion corresponds to the A = 4 operator F A F in the field theory, and will not play 
an essential role in our discussion below; 

(2) 

- the NSNS vector field B\ corresponds holographically to an operator of dimension- 
ality A = 2 + a/9 + 4Qj = 5 + 0(Qf) and sits in the supermultiplet with the term 
Tr(W d W^) + ---; 

- the fields C\ and A\ correspond to a transformed basis for two vector fields holograph- 
ically dual to Aj R = 3 and A = 7 operators. The operator Jr generates the U(1)r 
symmetry [29j [31] ; 

- the operator dual to the scalar field .Ao is sensitive to the value of a. For a = 1 (anti- 
D7-brane case) the scalar is dual to an operator of dimension A = 6 + Qf. For a = — 1 
(D7-brane case) the scalar is dual to an operator of dimension A = 2 + Qf [28]. In the 
probe approximation its mass sits on the AdS^ BF bound and the scalar is dual to [9J 

O = q ]a a a ^ , (5.5) 

where q is the doublet of SU(2)r squarks and a a p are Pauli matrices. The scalar 
operator is a vector of SU(2)r, which is dubbed "R-spin" in [6j. The case a = — 1 is 
the one we will focus on. 

In appendix [2] we have redefined the gravity fields dual to the A = 5, 7 operators into 
gauge-invariant combinations (which we dubbed t]b and r\c respectively, but see the comments 
around (|A.17|) and (|A.27[) ). which allowed us to decouple completely the e° equations. The 
scalar dual to the A = 2 operator (dubbed t\a) is sourced by the former gauge invariant 
combinations (as in equation (|5.2p ). 

Mixing of perturbations' equations of motion is a generic situation in holographic systems. 
In [3D] this situation was studied in generality and it was found that the generalization of the 
two-point function's prescription of [UJ consisted, schematically, on the expression 

G R oc B ■ A' 1 , (5.6) 
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the fields cj 0) 


= X, VA, 
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*x+-<pa Vx^wb ^x^-vc ' 
Vm-<pa Va^pb Va<-<pc 
Vb<-<p a V b <- 9b Vb^ Vc 
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. 
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where the subindex # indicates that we are talking about the retarded two-point function 
(which is imposed by fixing incoming wave boundary conditions at the horizon |41j). B is 
schematically a matrix of normalizable modes (dual to vevs of the operators [3D]) and A a 
diagonal matrix of non- normalizable modes (dual to sources of the operators). Focusing on 



(5.7) 



where (px is the source of the operator dual to the field X, and Vx^-^y is proportional to the 
vev of the operator dual to the field X when only <py is non-zero. This, together with the 
ingoing-wave condition at the horizon, determines the boundary conditions used to calculate 
the Vx^ VY , which are presented in the next section. 

The position of the quasinormal modes is a property of the retarded two-point function 
as a matrix, as the poles in a meromorphic expansion, and therefore all the Vx<^<p Y terms, 
are described by the same set of QNMs. 

In the former discussion we have not considered possible divergences arising in f)5.6|) . i.e., 
we have assumed that a holographic renormalization procedure has been made such that we 
can identify the matrix B in that expression with the matrix of normalizable modes. We have 
not proved that a proper cancellation occurs when explicit counterterms are added, especially 
with the lack of an asymptotically locally AdS spacetime. Here we assume this is the case. 
Indeed, background subtraction is enough to obtain the correct thermodynamical relations 
in which the free energy is obtained via the on-shell action, which is UV divergent. That 
the normalizable mode is the relevant factor to calculate the two-point function is a well 
established fact in holography [42] , even in the presence of irrelevant operators [43J o 

5.2 The retarded Green's function: numerical results 

The fluctuation of the axion at leading order (|A.ip was already studied in the literature 
J, and the quasinormal modes were found to be at w = 27rTn(±l — i) with n > an 



integer. In particular the frequency- function x ^ Vx can be found in [H]. As we will see below, 

the rest of the components — t^£x_ we ana ly Ze d no t show QNMs at oj = 27rTn(±l — i), even 
when the QNMs are shared by all components of the Green's function. The reason for this 
is the e* expansion. If we had included higher order corrections, the couplings between the 
equations of motion would have given rise to the presence of QNMs at precisely these values 
(plus an e* correction). In different words, the residue of these poles for the other — t^fx. 
components is subleading in e*. 

The only physical parameter appearing in the leading order equation of motion for the 
fluctuation rjc, (jA.llj) . is the temperature, which is encoded in the radius of the horizon 



21 Due to the structure of the effective perturbation action, counterterms contribute only to the hermitian 
part of the Gr matrix. Therefore, for the two-point function the anti-hermitian combination 

p(u) = i (Gr(u) - G«(w) + ) , (5.8) 

should not be sensitive to the addition of such counterterms. In [3D] it was shown that the spectral function 
matrix can be defined at any radius, contrary to the two-point function which is defined at the boundary. 
Since by construction all fields are regular at the horizon, the definition of the spectral function at the horizon 
proves its finiteness. The regularized hermitian part can be obtained by use of Kramers-Kronig relations. 
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Figure 4: Contour plots of the real (left) and imaginary (right) parts of 
oj/rh frequency plane. The real (imaginary) part is even (odd) under ui 



■pc 



ye 



in the complex 



-uj. 



rh = irT + 0{e*). This means that at leading order in e* we can work with the dimensionless 
ratio oj/rh- Taking r^ = 1 for convenience we present in figure H] a plot of the real and 



imaginary parts of 



V C+ 



fa 



<Pc 



in the negative imaginary frequency plane (no non-analyticities 



were observed in the upper half plane), following appendix [A] As explained in the appendix, 
the contour plot shown in figure 0] is valid up to some contact terms which, in particular, 
will not introduce any non-analyticities. Therefore we can say that at leading order in the 
backreaction parameter, for frequencies not much larger in norm than irT, there are no QNMs 
in the propagator of rjc- 

Similar comments hold for — ^7 VB , which we present in figure up to contact terms, 
for the negative imaginary frequency plane (no non-analyticities were observed in the upper 
half plane). We see now the existence of a quasinormal mode in the imaginary axis for 
uj/rh ~ — 2. hi. This mode is purely damped and gapped, therefore it does not correspond to 
a hydrodynamic mode. The fact that this QNM is not observed in the perturbation for the 
axion is, once again, an effect of the e* expansion. On the other hand, we will see below how 
this mode appears in other components of the Green's function. 



Now we come to 



V A^ VA 
9A 



. As seen in (JA.30J) there are two dimensionful physical parameters: 
the temperature (related to ru) and the charge density (related to rd). As we argued in the 
previous section, our solution is valid only for temperatures larger than the charge density 
per quark, so it is natural to scale dimensionful quantities with the temperature, namely oj/rh 

and r d /r h ~ n g ' jT . 



In figure we show a typical contour plot for 



V/u 



PA 



VA 



where three QNMs can be observed. 



One of the QNMs corresponds to a purely damped mode and presents a gap (its imaginary 
part is non-zero). In figure [7] we show how the purely damped mode behaves as we change the 
charge density relative to the temperature. In this plot we show the curve (r^/r/i) -3 as a guide 
to the eye, however we find a remarkable agreement between this line and the actual position 
of the damped mode, although we did not find an analytic explanation for it. Assuming this 
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Figure 5: Contour plots of the real (left) and imaginary (right) parts of 
oo/rh frequency plane. The real (imaginary) part is even (odd) under ui 



V B^ VB 



■rB 



in the complex 



-uj. 



mode is actually located at 



n 



(5.9) 



would imply that the mode will never become unstable by crossing to the upper-half plane, 
it will just approach the real axis as the radius of the horizon is small respect to the charge 
density. This is what was found in [Bj , where a zero sound mode is found at the origin for the 
rh = case (T = in the probe limit). 

The behavior of the oscillatory modes (those with a real part in the complex frequency 
plane) with the charge density is given in figure [8j where it is observed that increasing the 
ratio rd/rh pushes the modes deeper into the complex frequency plane (where the frequency 
is given in units of the radius of the horizon) . 

Sourcing explicitly the operators dual to t]b and r\c sources a vev for the operator dual to 
77,4, as can be seen from equation ()A.30p . where the presence of non-trivial ijb, rjc would trigger 
a non-trivial profile for A, even if we keep the boundary condition (pj{ = 0. We indeed observe 
this when we integrate the equations of motion with boundary conditions (a) ip^ = <pb = 
and (pc = 1 or (b) 994 = <pc = and ps = 1- hi the case (a) the contour plot for — ^~ VQ is 
given by figure [H and the one for — ^ Q in figure [9] (left). We can see quasinormal modes 



pc 



whose positions depend on r^jr^. The position of these quasinormal modes is the same as in 

:ZB - (see 



figure El In the case (b) again — ^* B is given by figure El and the contour plot for 



figure [9] (right)) presents quasinormal modes with positions coincident with the quasinormal 
modes observed both in figure El and HI These coincident QNMs are not surprising, as we 
commented already, since these equations are coupled and QNMs are shared by the different 
components of the matrix of Green's functions. This is the dual to the mixing of operators in 
the field theory side, such that sourcing one single operator inevitably implies a vev for more 
than one operator. 

Summarizing the results presented in this section, we observe that all the QNMs at leading 
order in e* are in the lower complex plane, and therefore are stable modes. Going to next 
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Figure 6: Contour plots of the real (left) and imaginary (right) parts of 



<PA 



in the complex 



oo/vh frequency plane for r^jvh = 1.4. The real (imaginary) part is even (odd) under u 
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Figure 7: Position of the purely damped mode in the oo/r^ imaginary axis. The dashed line 
corresponds to (r^/r^) -3 . 



order in e* would require to solve equations of motions where the corrections to the back- 
ground solution found in section [2] appear explicitly, however such a correction could modify 
the position of the quasinormal modes only perturbatively in e*, and would not trigger an 
instability. The only case where such a correction might have an important effect would be in 
the rhjr<i — > case, where (|5.9p points to a mode sitting at the origin of frequencies. However, 
this situation is out of the regime of validity of our solution. 

6 Conclusions 

We have expanded the study of the charged D3-D7 black hole started in [24J providing a new 
solution for generic value of the parameter r<j. The solution is dual to M = 4 SYM (and quiver 
generalizations) with massless dynamical flavors at finite temperature and baryon density. Let 
us discuss in turn two features of this solution. 
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Figure 8: Position in the oj/r^ complex plane of the first quasinormal mode with non-null 
real part (by symmetry there is another QNM related to the one showed in the graph with 
opposite real part) from r^jr^ = 0.6 (point A in the graph) to r^jr^ = 1.6 (point B). At lower 
values of r^/r^ < 0.6 the position of the quasinormal mode does not change perceptively. 



6.1 Breaking of the perturbative solution and of the probe approximation 

From equation (|3.T|) (using n q = N c ub and specializing to the flavored N = 4 SYM case), we 
observe that the dimensionless quotient r^/r/j is a measure of the charge (baryon) density per 



flavor, n B / q = jf-, over temperature 



r|_ *n B/q 

ri ~ va:t3 + U{e *> ■ [bA) 

From this formula it follows that, since r^/r^ can be of order one in our solution, the latter 
can describe effects due to large charge density per flavor (order y/X^T 3 ). 

Nevertheless, the overall value of the charge density is still "small", in the sense that it 
is not sufficient to modify at order one the IR physical observables. In fact, the latter are of 
the generic form O ~ Oj\f=i [1 + e/(r^/r/j) + ...]. As we have seen in section [H e.g. for the 
temperature, the flavor corrections e/(r<j/r/i) are of order one when (rd/r^) 3 ~ 1/e, but, as 
expected, this value makes higher order corrections in e relevant, thus spoiling the reliability 
of the solution. 

From the gravitational side this is expected since the gravitational force would not be able 
to compensate the electric repulsion in a small black hole. It is not clear to us what to expect 
in the field theory side. A transition to a different phase, like the spatially modulated phases 
encountered for the dual of the AiS-RN black hole [47J would appear if instabilities are found 
at finite momentum, which we have not covered in the present work. 

Thus, while the solution can correspond to large charge with respect to the temperature, 
it does not describe in the IR an (almost) zero temperature system, which would require a 
solution at all orders in e/^/r^). The existence and the nature of such solution remains an 
extremely interesting open problem. 

The breakdown of the perturbative solution at zero temperature implies the breakdown 
of the charged probe approximation as well. In fact, the probe approximation requires that 
the original background is not significantly affected (or it is un-affected for a long time) by 
the addition of a small perturbation due to flavor branes (i.e. tf{rd/ r h) "C 1). But in the 
present case the original background does not admit a small perturbation at zero temperature: 
either the perturbed system is at non-zero temperature (and "far" from zero temperature in 
a qualitative way, as we have seen in this paper), or the flavor contribution must be dominant 
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case (b) discussed in the text, with X = B. In these figures ra/vh = 1. 



near the origin of AdS in order to have a zero or almost zero temperature (the corresponding 
solution being still unknown). 

This breakdown of the charged probe approximation is actually already known. For exam- 
ple, in [5] it is pointed out that for the probe approximation to make sense, the stress-energy 
tensor of the probe must be subleading with respect to the one of the background, otherwise 
one has to take into account the brane backreaction. In a zero temperature background and 
for a charged probe, this condition cannot be satisfied all the way up to the deep IR region! 22 ! 
unless there is some other mass scale, such as the temperature or the flavor mass. In the latter 
cases, the probe approximation is valid as long as the charge density is not parametrically 
larger (in N c ) than these extra scales. 

6.2 Stability 

With the aim of finding a physical instability we have checked the matrix of susceptibilities 
of the system, which turned out to be parametrically stable, and quasinormal modes of a 
set of perturbations of the system, showing no indication of a physical instability either. We 
comment further on the quasinormal modes below. 



Comparison of the fluctuation calculation with the probe limit. In equations (J5JJ)- 
(|5.3|) we have sketched how the equations for the fluctuations organize when a perturbative 
expansion in e* is implemented. In the probe approximation fluctuations of the supergravity 
fields are not allowed and Ao = 0. This satisfies identically equation (|5.ip and transforms 
equation (|5.2p into the probe fluctuation equation EOM[X\\ = 0. A non trivial solution 
to these equations will source fluctuations of the supergravity fields Ai via equation (|5.3p . 
which are first order in e*, which is a depreciable correction to the background value in the 
probe limit. It is only after a large time of order t ~ to/ e * (*o being order one) [H] that the 



22 Condition (3.27) in [5] is equivalent to the condition for the validity of our solution (rd/rh) -C 1/e. 
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correction Ai would not be small compared to the background value and the backreaction of 
the probe onto the geometry needs to be considered. 

So we see how the e* expansion of the fluctuations organizes itself to take into account the 
backreaction of the geometry, and how can one link the results to the probe approximation 
case. The probe calculation ceases to be valid at t ~ io/e* where perturbations of the probe 
on the horizon must be considered: energy is injected into the system and eventually it will 
be absorbed by the black hole, which will increase in size. 

In the specific case of the fluctuation of the field .4.0 ( see section IA.4J) , the fact that t]b 
and r)c source the equation of motion for t]_a points to the fact that the source will behave 
differently when embedded in the original background or in the perturbed one. However, if 
we turn off the source of the irrelevant operators that perturb the background, the probe will 
be insensitive to these changes (to leading order in e*), and we obtain the same equation of 
motion one would have obtained from the probe calculation at finite temperature. 

Operator mixing and the two-point function. The equations of motion for the fluctua- 
tions we have focused on consist of a set of coupled, linear, second order differential equations. 
In section [5] we have expanded these equations in a power series in e* , since our background 
is only valid up to first order in that expansion. This allowed us to decouple the e* equations 
for the fields dual to the A = 5, 7 operators, however the field dual to the A = 2 operator, 
itself of order e*, is sourced by the former ones. 



Equation (|5.6p describes the mixing of the operators in the Green's function. We have 
given the values of — < ~ y,A , — t^M. anc j — tz^o. j n figures [UJO and [3J respectively, and x ^ ipx can 

be found in [13]. The off-diagonal terms, which describe the mixing of operators, — t^£B_ anc [ 

— — ^ correspond to the values in figure We have observed that the position of the QNMs 
in the off-diagonal two-point components coincide with the positions of the QNMs in the 
two correspondent diagonal components, as expected from operator mixing considerations. 
We have not described all of the Vx-tr-ip Y components, since we were mainly interested in the 
gapped, purely damped mode described in figure We do not expect that the inclusion of 
the rest of the components will add any relevant physical information. The main interest of 
considering all the components would be to check the relations that are imposed on the two 
point matrix function by time reversal symmetry 

G R (-u) = G R (u>)* , G A (u) = G R (u,y . (6.2) 

Absence of overdamped modes. The fluctuation of the backreacted system shows no 
sign of instability at zero momentum. At finite temperature the QNMs have a negative 
imaginary part and e* corrections, being parametrically small, cannot push these modes to 
the positive half-plane. The reason for this is that our solution, as stated several times now, 
holds only when the temperature dominates over the charge of the system, or when the two 
effects are of the same order. 

We have focused in a purely damped mode with a gap at finite temperature apparently 
given by equation (|5.9j) . In the probe case at zero temperature this mode has been observed 
in [6]. For this mode to be present in the probe calculation it is crucial that fluctuations of 
worldvolume fields effectively see an AdS2 spacetime in the IR, which is associated to extremal 
black holes (and therefore zero temperature solutions). In the present work we expect the 
zero temperature setup to be described not by r^ = 0, but an extremal black hole with finite 
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horizon radius, whose near-horizon geometry would be AdS2 X R 3 , and this would imply that 
the gap is non vanishing at zero temperature. 

We have shown that the low temperature regime is beyond the applicability of our solution, 
so the former comments are just speculative. Even if we had a T = solution at hand it 
is possible that no instability would appear. The reason is that in [28] it was shown that 
backreaction corrects the mass of the field dual to the A = 2 operator to a larger value, getting 
away from the BF bound in the stable direction, the dual operator having now dimension 
A = 2 + Qf > 2. Therefore, we can speculate that if instabilities of the system were to be 
found they arise as new solutions, presenting a spatially modulated phase or maybe something 
more complicated. 
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A Expressions involved in the study of the perturbations 

By construction, the equations of the modes we focus on are those that vanish identically 
in our background. Specifically, these are the equations of motion for the field strengths 
F± , Fj , Fq j, F 7 1 , H I , fig and dA\ , which are a total of twelve equations so that our 
system may seem overdetermined (see |28j). However, it is not difficult to show that the 

(3) 

equation of motion for F$ still vanishes identically, reducing the number of equations to 
nine and suggesting that the system is actually under-determined. This conclusion is also 
too naive. The reason is explained in [28], where it is shown that acting with the external 
derivative on the equations of motion for F^ and H2 gives rise to the equations of motion 
for F^ and dAo, respectively. The ultimate reason for this is that the modes Cq and Bq 
are nothing but Stiickelberg scalars coupled to the vectors C[ and B\ o The equation of 

(5) 
motion for F^ does not provide any new physical information, since it is not independent of 

the other equations of motion. On the other hand, the equation for dAo does provide a new, 

independent equation. 

A.l Dimension 4 scalar operator 

The first equation is the one for the fluctuation of the axion and at leading order 0(e*) is 
given by the equation of motion of a massless scalar in Schwarzschild-^4dS* 

cf + d r log [r(r 4 - ri)] cf ' + j~^C^ = . (A.l) 

v' 'h> 

23 Actually, the Stiickelberg scalar associated to the vector B[ is a combination of the scalars B and Ao- 
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A. 2 Dimension 7 and 3 vector operators 

The second set consists of five different equations for the fluctuations of A\, C\ and Cq . At 
leading order in e* these are fluctuations of N = 4 SYM that can be considered independently 
of the backreaction of fundamental matter. It is convenient to express the vector fields in 
a different basis where we have a massive and a massless vector field dual to operators of 
dimension As = 7 and Aa = 3 respectively, which is achieved by taking 

__ -y / 2A + 2E ( 4 )_ 2A + v / 2S 

Ai- V2 , C x - - . (A.2) 

Then the five equations for the fluctuations read 

w(wA r - iA' t ) = , (A.3) 

3wA r + ru;A' r - i(3A' t + rA'{) = , (A.4) 

-3i(r 4 - rfl)uj^ r + 24r 3 S t + 24ir 3 wC^ 4) 

-i(r 4 - rt)(rcoK - i(3^ + rE?)) = , (A.5) 

(5r 8 - 6r 4 r£ + r£)£ r + ir 5 wS t + r(r 4 - r A h ) 2 T,' r - r 5 w 2 C^ 4) 

_ (5r 8 _ 6r 4 r 4 + r 8) C W' _ r(r 4 _ r 4 )2(7 W" = Q ^ (A _ 6) 

-\/2r 2 cj 2 A r + (-48r 4 + 48rj| + 2r 2 w 2 )£ r + ir 2 u(V2A' t - 2Ej) 



+48(r 4 -r 4 )ci 4) ' = 0. (A.7) 

Notice also that At enters only via derivatives in the equations of motion, so we can always 
make a shift to set the constant part of At to zero. 
To solve (JA.3P and (|A.4j) we just need to impose 

A r = -A' t , (A.8) 

U) 

where we are assuming cj ^ in the following. With this plugged back into the equations, all 
dependence on the At r functions disappears, and one can algebraically solve for S r as well 



j)cf - *r 2 
24(r 4 — rt) — r 2 uj 2 



24(r 4 - r 4 )d 4) - ir 2 ujZ' t /A . 

^ = \.^ h ! _ 2 .. 2 . (A.9) 



and plugging ()A.9|) into (|A.5|) and ()A.6|) one gets the same differential equation involving 
zeroth, first and second derivatives of St and Cq . The solution to this equation is given by 

Cf = -£t + ^, (A.IO) 



where rjc satisfies 

rjlj + d r log 



r 3 (r 4 — rf) 



24( r 4 _ r 4) _ r2a; 2 



24r 2 (r 4 -r 4 )-r 4 w 2 



Vc { — 4 % VC = • (All) 



Before solving this equation of motion let us outline that the relations (|A.8j) and (JA.9P are 
ambiguous. The fact that we cannot fix completely the five modes present in these equations of 
motion is reflecting the existence of a gauge freedom. In particular, two modes can be gauged 
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away (this can be seen from the set (|A,3p - (|A.7j) by noticing that there are two differential 
relations between these equations). 

To solve numerically the equation of motion for the fluctuation it is better to work with a 
compact radial variable defined as z = 1/r. In these coordinate I 24 ! equation (|A.lip is solved 
near the boundary z = by 



s c s c u} 2 s c ^J 4 s c uj 2 4 2 scw 4 (1728r£ + 



u 



^ "^ + ^4^ + TToY " r *"° " IT V 442368 Z + Vc * (A ' 12) 

w 2 (-3317760v c + s c w 2 (262656r 8 + 1848r 4 w 4 + w 8 )) 8 
+ 53084160 Z + '" 

/ ^ 2 (147456r 8 t + r728r 4 V + ^ 6 
+ ° g[ J V 2211840 

s c ^ 4 (147456r 8 + 1728r 4 w 4 + w 8 ) 8 

~Z i 



35389440 

and we see that the asymptotic solutions in the UV behave as z~ 4 and z 6 , corresponding to 
the dual to a vector operator of dimension A = 7. 

In |43j it was shown that for irrelevant operators the coefficient multiplying the normal- 
izable mode, proportional to z 6 , is still holographically dual to the expectation value of the 
operator, whereas the non-normalizable z~ 4 mode is dual to the source of the operator in the 
field theory. 

Near the IR there are two independent solutions, one corresponding to a ingoing wave at 
the horizon 25 ! and one corresponding to an outgoing wave. Focusing in the ingoing wave we 
obtain 



[ 8u(2ir h + oj) \ r h 

where f]c,h is an undetermined normalization. 



(A.13) 



We consider now regularity at the horizon, which is obvious from the previous expression 
for rjc, but not evident for the physical modes Cq and Sj. A quick look at the previous 
equations shows that, assuming regularity for the time components, the only problems that 
can arise come from equation (|A.9p , Using an ingoing wave regular solution for the time 
component of the gauge field, and the r\c expansion ()A.13p . shows that to have a regular 
radial component of the vector field, Si r , the fluctuation of the time components must vanish 
at the horizon, i.e., 

^t = #{r-r h f-^ +••• . (A.14) 

Imagine for a moment we pick the Stiickelberg gauge in which Cq = 0. Then from (jA.lOp 
we obtain that the time component Sj has a finite value at the horizon (proportional to rjc h) 
and the radial component would be irregular at the horizon, behaving like S r ~ (r — r/j) -1 . 
This is therefore an irregular gauge. From now on we will pick the regular, radial gauge 26 ! 



24 In the rest of the appendix we will use z radial coordinate in expressions necessary for the numeric 



integration and r radial coordinate otherwise. 

25-t 

26 The same discussion will hold for the fluctuation of the NSNS potential below. 



25 We are assuming here that the temperature is non-zero, rh 7^ 0. 
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S r = 0. We obtaiitJ 

f , 24i(r 4 - ri)uj 
S *= ldr n4f J ^ h) ^J r (A.15) 



24( r 4 _ r 4) _ r 2 w 2 ' 



and performing an integral we obtain the complete function. When one does this integral 
a constant appears that must be fixed somehow. This constant is actually unphysical. In 
the equations of motion the terms proportional to time components of the massive vector 
fields without a derivative appear always in a certain combination with the corresponding 
Stiickelberg scalar without derivative, and the constant shift of the time components can 
be compensated by the opposite shift in the Stiickelberg scalar. We can use this to set the 
constant of integration in the previous equations to zero. 

The only physical parameter appearing in the equation of motion for the fluctuation rjc 
is the temperature, which is encoded in the radius of the horizon r^ = ttT + 0(e*). This 
means that at leading order in e* one can build the dimensionless ratio co/rh and describe 
all the finite temperature solutions at the same time, a change of temperature scale being 
compensated by a change in frequency. 

From equation (|5.6p we are interested in determining Vc-^-w c and <pc- Prom the holo- 
graphic correspondence the last identification corresponds just to the boundary condition 
Vc = sc, whereas the former must read Vc*-sp c oc vc [13] when sc is non-zero. We are not 
interested in determining the exact proportionality factor, but in the position of the poles. 
Therefore, to produce the plots in figure 2] we solve equation (|A,11|) performing a double 
shooting from the IR and the UV with sc = 1 in the boundary conditions (|A.12|) and (|A.13p , 
and varying 7]c,h an d vc such that the solution and its derivative match smoothly in the 
bulkH Then 

^^ oc ^ , (A.16) 

Vc s c 

and we will renormalize the Green's function such that the constant of proportionality is 
exactly one. Actually, to study the two-point function we are interested in we should not 
focus just on vc/sc but in the result coming from undoing the change of variables (|A.15|) . 
Since we are evaluating an antiderivative (no extra constant must be considered) it can be 
performed with the asymptotic series in the UV, to obtain 

(0^)n oc *L _!>!_!>! _ " 10 . (A.17) 

N ' s c 36 3072 5308416 v ' 

The frequency powers correspond to contact terms (derivatives of Dirac's delta distributions 
when we Fourier-transform back). In particular they are analytic in the complex frequency 
plane. Our main focus is to determine the non-analytic behavior of the two-point function, 
so we may ignore these contact terms. 

A. 3 Dimension 5 vector operator 

The third set of equations of motion contains three coupled equations. As in the previous 
subsection, at leading order in e* these are fluctuations of M = 4 SYM that can be considered 



27 Notice that the solution for A, i.e., the field dual to the Aa = 3 operator, will not be needed; actually, 
this solution cannot be found uniquely since this is a pure gauge mode in the case of zero momentum. 

28 The numeric condition we have set for this smooth juncture is that the maximum difference between the 
function or its derivative at the matching point is less than 10 -10 . 
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independently of the backreaction of fundamental matter 

8r(r 4 - ri)B { $ - 2^2r(r 4 - r£)<$sj 2) ' - r 3 u(ujB^} - iB$) = 
2V2ir 3 iuBl 2} - 3i(r 4 - rf^coB^ + 8r 3 s{ 2) - iur(r A - rfi)B$' 

-3(r 4 - r 4 )^' - r(r 4 - r 4 )BJf = 
^J 2) + 2V2ir 5 uB ( y ) + (5r 8 - 6r 4 r£ + r 8 h )(2V2B {2) - B {2) ') 

+r(r s - 2r 4 r 4 + r 8 h )(2V2B[ 2 / - B {2) ") = 



-rV« (2) 



As in the previous subsection, there is an algebraic solution for the radial component 



B 



(2) 2^2(r 4 - rj)Bf - ifojB^ 



lr 



8(r 4 — rf) — r 2 cj 



2,, ,2 



(A.18) 
(A.19) 
(A.20) 

(A.21) 



which, once plugged back into the equations of motion, identifies (JA.19P and (|A.20j) in an 
equation with solution 



(2) 



with rjB satisfying 

T\ B + 9 r log 



B, 



np [ , y a if** 



2^2i 



to 



B[?+VB, 



8(r 4 — rf) — r 2 uj 



8r 2 (r 4 — r 4 ) — r A uj 2 

^ (r 4 _ r 4 )2 VB = 



(A.22) 



(A.23) 



As before, there is an ambiguity due to the gauge freedom (which can be seen in the set (JA.18J) - 
(IA.19P by noticing the existence of a differential relation between the equations). Equation 
(JA.23P has as asymptotic solution in the UV the expansion 



sb , s b uj" 2 i / v B u 2 t 7s B cj 4 (192r^ +W 4 ; v ,. : 
VB W-J H 7S r ft S S Z + V B Z + I 1 TT^T I z° - 



16 



+ log [z] 



s B w 2 (102r 4 +w 4 ) 4 s B w 4 (192r 4 +w 4 ) ( , 



768 



-** + 



6144 



49152 

z 6 + • • • 



(A.24) 



where the term proportional to the non-normalizable mode, sb, is dual to the source and the 
term proportional to the normalizable mode, vb, is related to the vev of the dual, dimension 
A = 5 operator. Near the IR the ingoing wave solution is given by 



7] B ~ (1 -r h z) 4r hrj By 



128r s h - 32ir 2 h cj + 6r h oj 2 + 2,iuj A ( 1 
8u}(2ir h +u) V r h 



+ ■■■ 



(A.25) 



with r]B,h a normalization. 

The same discussion about regularity at the horizon that was performed for the £ vector 

(2) 

field can be performed in the present case, and we pick the radial gauge Br = to ensure 
regularity. Then we have the relation 



B 



(2) 
It 



dr 



2y/2i( 



r — rZ)u) 



tVb 



8(r 4 — rf) — r 2 co 2 
where the integral must be considered, once again, as an antiderivative. 



(A.26) 
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(A.27) 



Once again, the only physical scale at leading order in flavor backreaction is the temper- 
ature, and we can express the solution in terms of uj/r^. The results presented in figure [5] 
correspond to the identification 

Vb<-ipb __ vb_ 

PB S B 

which can be obtained with a double shooting from the IR and UV with sb = 1 and the 
parameters r]B,h an d vb adjusted to match smoothly the solution and its derivative. The real 
two-point function component we are interested in must be obtained undoing the transfor- 
mation (|A.26p . which gives us a result 



(O b O b )r<x — 

SB 



rfa 2 



CO 



1024 ' 



(A.28) 



and again the frequency powers correspond just to some contact terms that we can ignore 
since we want to focus on non-analytic terms. 

A. 4 Dimension 2 scalar operator 

The modes studied in the previous subsections enter as sources in the first order equation of 
motion for the fluctuation of the worldvolume scalar, <L4o- Plugging the previous algebraic 
solutions we find that the equation of motion for the worldvolume scalar is solved by 



with r/,4 satisfying 



A 



-R (2) +r7„ 



(A.29) 
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2r 10 + 2r 6 r 4 
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rj' A + d r log 



/y*\J I fy*^ 



4r 4 



r h) + ( r6 + r , 



6 W 
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>U 



VA 



r 4 r 6 d + hr\ri 

Vb 



r ( r 6 + r 6)( 8 ( r 4 _ r j\ _ r 2 w 2) 



V2- 



(r 6 + 2r\ + 2r 3 



r b + r c 



-2i 



(r 6 + r 6 d )(r 4 - r 4 ) 
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(r 6 + r6)(24(r 4 -r 4 ) - r 2 u 2 ) 



-VB 



V'c 







(A.30) 



Since Aq is dual to a dimension A = 2 operator the two independent solutions in the UV 
behave as z 2 and z 2 log[z]. In these limiting cases where the Brietenlohner-Freedman bound 
is saturated, it is the coefficient of the most divergent term, namely the one with z 2 log[z], 
the dual to the source of the operator, whereas the coefficient of z 2 will be proportional to 
the expectation value [16]. Since t\b and r\c are dual to irrelevant operators that deform the 
UV of the theory, the asymptotic behavior is dominated by these fields [33]. Using (|A.12p 
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and (|A.24|) . we find 



VA 



SB 



2V22 



uj(Wiscr d — 3\/2s_b^) 9 

192 * 



(A.31) 
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55ir 3 d s c oj 5 + 288(4\/2vb + (-4s^ + y/2s B r A h - <iv A )uj' 

4608 



-, 4 + 




2 a; 2 (-768s^ + 48is c ?> 3 + V2s B (192rf i + w 4 )) , 



3072 



z^ + 
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128 



(A.32) 



The presence of the logarithm introduces a scale zo (which we fixed to 1) that enters in the 
counterterms one has to add to regularize the on-shell action. This affects the definition of 
the two-point function introducing contact terms (analytic terms in the momentum), which 
we are ignoring because we are interested only in the position of the poles of the retarded 
two-point functions, so our results will be insensitive to this scale. For the IR, the solution 
with the ingoing wave boundary conditions is given asymptotically by 
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To integrate the equation, we perform again a double shooting from the IR and UV. Since 
we are interested in the two-point function associated to the operator of dimension A ^4 = 2 
we must make sure that we are sourcing only this operator. This implies that we have to 
integrate with conditions s.4 = 1 and sb = sc = 0. From this integration we can extract 



V A < 



-<PA 



VA 



VA 
SA 



(A.34) 



The result for va/sa is reported in figure 



The operators dual to rjs and rjc will source as well the operator dual to tja- This is seen 
from equation ()A.3ip , where the presence of non-trivial t]b or r/c would trigger a non-trivial 
value of va, even if we keep the boundary condition s A = 0. We indeed observe this when 
we integrate the equations of motion with boundary conditions (a) sa = sb = and sc = 1, 
giving 

Vm-jpo __ VA 
>~Pc s c 

or (b) sa = sc = and sb = 1, giving 
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